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. Abstract 

(N 

A geometric generalization of discrete-time linear deadbeat observer 
is presented. The proposed method to generate a deadbeat observer for 
a given nonlinear system is constructive and makes use of sets that can 
. be computed iteratively. For demonstration, derivations of observer dy- 

namics are provided for various example systems. Based on the method, 
a simple algorithm that computes the deadbeat gain for a linear system 



■ with scalar output is given. 



> 



1 Introduction 



Observer design for linear systems is generally acknowledged to be understood 
well enough. For discrete-time linear system — Ax with output y — Cx, 
, Luenberger observer [5] dynamics read 

^ ; x+ = Ax + L{y - Cx) (1) 

o ■ 

' and designing the observer is nothing but choosing an observer gain L that 

places the eigenvalues of matrix A — LC within the unit circle. Simple and 
elegant, anyone would hardly doubt that this construction is the construction 
for linear systems. However, perhaps due arguably to over-elegance of the nota- 
' tion, it is nontrivial to unearth the true mechanism (if it exists) running behind 

Luenberger observer in order to generalize it in some natural way for nonlin- 
ear systems. In this paper we aim to provide a geometric interpretation of the 
righthand side of ([T]) for the particular case where matrix A — LC is nilpotent, 
i.e., when the observer is deadbeat. Our interpretation allows one to construct 
deadbeat observers for nonlinear systems provided that certain conditions (As- 
sumption [T] and Assumption hold. We now note and later demonstrate that 
when the system is linear those assumptions are minimal for a deadbeat ob- 
server to exist. The literature on observers accommodates significant results. 
See, for instance, [5l|4l|9l[l3l[IIl|3l[Il[l4]. 
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The toy example that we keep in the back of our mind while we attempt to 
reach a generalization is the simple case where A is a. rotation matrix in 



A = 



cost 
sin^ 



- sin 61 
cos 6 



with angle of rotation 9 different from and tt. Letting y — X2, i.e., C = [0 1], 
the deadbeat observer turns out to be 



X 



Ax - 



cos 29/ sin 9 
sin 26*/ sin 6* 



(y - Cx) 



which can be rewritten as 



cot 6* 
1 



[y - Cx) 



Now we state the key observation in this paper: The term in brackets is the 
intersection of two equivalence classes (sometimes called congruence classes [B]). 
Namely, 



cot 6* 
1 



{y - Cx) ^ {x + Anull(C)) n (a; + null(C)) 

as shown in Fig. [1] Based on this observation, one contribution of this paper 

,i + Anull(C) 



x + 




X + null(C) 



Figure 1: Intersection of two equivalence classes. 



is intended to be in showing that such equivalence classes can be defined even 
for nonlinear systems of arbitrary order, which in turn allows one to construct 
deadbeat observers. There is another possible contribution that is of more prac- 
tical nature: We present a simple algorithm that computes, for linear systems 
with scalar output, deadbeat gain L by iteratively intersecting linear subspaces. 
(Devising reliable numerical techniques to compute deadbeat gain for discrete- 
time linear systems had once been an active field of research; see, for instance. 
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The remainder of the paper is organized as fohows. Next section contains 
some prehminary material In Section [3] we give the formal problem definition. 
Section H] is where we describe the sets that we use in construction of the dead- 
beat observer. We state and prove the main result in Section [51 An extension 
of the main result where we consider the case with input = f{x, u)) is 
in Section ^ We provide examples in Section [71 where we construct deadbeat 
observers for two different third order systems. In Section [8] we present an al- 
gorithm to compute the deadbeat observer gain for a linear system with scalar 
output. 

2 Preliminaries 

Identity matrix is denoted by /. Null space and range space of a matrix M S 
T^mxn g^j.g denoted by J\f{M) and 7^(M), respectively. Given map ^ : X ^ y, 
yU~^(-) denotes the inverse map in the general sense that for y G y, H~^{y) is 
the set of all x G X satisfying = y. That is, we will not need /i be bijective 
when talking about its inverse. Note that y ^ A<('^) will imply fJ.~^{y) = 0. 
Linear maps x i— Mx will not be exempt from this notation. The reader should 
not think that M is a nonsingular matrix when we write M^^y. (In our case 
M need even not be square.) For instance, for M = [0 0] we have M~^y — 
for 2/ 7^ and M~^0 = R^. The set of nonnegative integers is denoted by N and 
IR>o denotes the set of strictly positive real numbers. 

3 Problem definition 

Consider the following discrete-time system 

x+ = f{x) (2a) 
y = Kx) (2b) 

where x G X <Z M" is the state, x^ is the state at the next time instant, and 
y G y d M™ is the output or the measurement. The solution of system ([2]) at 
time k G N, starting at initial condition x G X is denoted by (/)(fc, x). Note that 
0(0, x) — X and (l>{k 4- 1, a;) = /(0(fc, x)) for all x and k. 
Now consider the following cascade system 

x+ = fix) (3a) 
x+ G g{x, h{x)) (3b) 

We denote a solution of subsystem ([Sb) byi/'(fc, x). We then have -0(0, x, x) = 
X and 'ip{k+ 1, i, x) G gi^Pik, x, x), h{<p{k, x))) for all x, x, and k. We now use 
([3]) to define deadbeat observer. 

Definition 1 Given g : X x y ^ X , system 

G g{i, y) 
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is said to be a deadbeat observer for system ^ if there exists p > 1 such that 
all solutions of system ^ satisfy 

ipik, X, x) = 4>{k^ x) 

for all X, X X and k > p. 

Definition 2 System ([2|) is said to be deadbeat observable if there exists a 
deadbeat observer for it. 

In this paper we present a procedure to construct a deadbeat observer for 
system ([2]) provided that certain conditions (Assumption [T] and Assumption [2]) 
hold. Our construction will make use of some sets, which we define in the next 
section. Before moving on into the next section, however, we choose to remind 
the reader of a standard fact regarding the observability of linear systems. Then 
we provide a Lemma [T] as a geometric equivalent of that well-known result. 
Lemma [1] will find use later when we attempt to interpret and display the 
generality of the assumptions we will have made. 

The following criterion, known as Popov-Belevitch-Hautus (PBH) test, is an 
elegant tool for checking (deadbeat) observability. 

Proposition 1 (PBH test) The linear system 

x^ ~ Ax (4a) 
y = Cx (4b) 

with A e R"^" and C G R™^" is deadbeat observable if and only if 

= n for all A 7^ (5) 

where X is a complex scalar. 

The below result is a geometric equivalent of PBH test. 

Lemma 1 Given A e M"^" and C G M™^", let subspace Sk o/R" be defined 
as Sk '■= ASk-i n So for k = 1, 2, . . . with Sq := J\f{C). Then system (|4]) is 
deadbeat observable if and only if 

Sn = {0} . (6) 

Proof. For simplicity we provide the demonstration for the case where each Sk 
is a subspace of C" (over field C) . The case Sk C R" is a little longer to prove 
yet it is true. 

We first show ([6|) Suppose ([5]) fails. That is, there exists an eigen- 

vector V G C" and a nonzero eigenvalue A S C such that Av = Xv and Cv = 0. 
Now suppose for some k we have v G Sk. Then, since v is an eigenvector with 
a nonzero eigenvalue, we can write v £ ASk. Observe that v € Sq for Cv — 0. 



rank 



A -XI 
C 
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As a result v G ASk nSa — Sk+i- By induction therefore we have v E Sk for all 
k, which means that ([6]) fails. 

Now we demonstrate the other direction ([5]) ©• We first claim that 
Sk+i C Sk for all k. We use induction to justify our claim. Suppose Sk+i C Sk 
for some k. Then we can write 

Sk+2 = ASk+1 n So 

c ASk n So 

= Sk+1 ■ 

Since Si C Sq our claim is valid. A trivial implication of our claim then follows: 
dim5fc+i < dimiSfc for all k. Let us now suppose ([5]) fails. That is, dim5„ > 1. 
Note that dim5o < n. Therefore dim5„ > 1 and dim 5^+1 < dim5fe imply the 
existence of some £ G {0, 1, . . . , n — 1} such that dimSi+i — dim5^ > 1. Since 
5^+1 C Se, both Se+i and Si having the same dimension implies 5^+1 — Se- 
Hence we obtained Se — ASi D So which allows us to write Se C ASe. Since 
dim Se > dim ASi we deduce that Si = ASe ■ Since dim Se > 1, equality ASe = 
Se implies that there exists an eigenvector v E Si and a nonzero eigenvalue 
A e C such that Av = Xv. Note also that Cv — because Se C So- Hence ([5]) 
fails. ■ 

Remark 1 It is clear from the proof that if ^ fails then dhnSk > 1 for all k. 

4 Sets 

In this section we define certain sets (more formally, equivalence classes) asso- 
ciated with system For x G X we define 

[x]o ■.^h'\h{x)). 

Note that when h{x) = Cx, where C G M™^", we have [a;]o x+N{C). We 
then let for fc = 0, 1, . . . 

Wfc+i := [x\t n [x\o 

where 

[x]t:^f{[f-\x)]k). 

Note that [x\l = when x ^ f{X) since then f^^{x) = 0. 

Remark 2 Note that [x]k+i C [x]k and [x]^^^^ C [x]'^ for all x and k. 

The following two assumptions will be invoked in our main theorem. In hope 
of making them appear somewhat meaningful and revealing their generality we 
provide the conditions that they would boil down to for linear systems. 
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Assumption 1 There exists p > 1 such that, for each x E X , set is 
either singleton or empty set. 

Assumption [T] is equivalent to deadbeat observability for linear systems. Be- 
low result formalizes this. 

Theorem 1 Linear system is deadbeat observable if and only if Assump- 
tion]^ holds. 

Proof. Let Sk for fc 0, 1, . . . be defined as in Lemma [1] Note then that 
[x]o = X + Sq. We claim that the following holds 

( x + Sk for a; e TZ{A'') 
^^^'^ ^ \ for a; ^ n{A'') ^' > 

for all k. We employ induction to establish our claim. Suppose ([7]) holds for 
some k. Then we can write 

[x]+ - A[A-^x]k 

= A[A-^xnniA'')]k. 

Note that A-^x n 7^(A'=) 7^ if and only if x e n{A''+'^). Since [a;]fe+i [x]+ 
[x]o, we deduce that [x]k+i = for x ^ n{A''+^). Otherwise if a; G n{A''+^) 
then there exists some rj G Tl{A'^) such that Arj = x. Using this 7/ we can 
construct the equality A^^x — r] + J^{A) and we can write 

Nfc+i = Nfe n [x]o 

= A[A-\]k n [x]q 

= A[f]+U{A)]kn[x]o 

= A{rj + {U{A)mZ{A''))+Sk)n[x]o 

= {At^ + ASk) n {x + So) 

= {x + ASk) n {x + So) 

= x + {ASknSo) 

= x + Sk+1 ■ 

Since ([7]) holds for k = 0, our claim is valid. 

Now suppose that the system is deadbeat observable. Then by (O we see 
that Assumption [1] holds with p = n + 1 thanks to Lemma [TJ If however the 
system is not deadbeat observable, then by Remark [T] dim 5^ > 1 for all k. We 
deduce by ([7]) therefore that [0]^ can never be singleton nor is it empty. Hence 
Assumption [T] must fail. I 

Assumption 2 Given x, x E X and k; x G [x\^ implies [x\l = Wi^ ■ 

Theorem 2 Assumption]^ comes for free for linear system (|3|). 
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Proof. Evident. ■ 

Last we let [x]'ti X and define map n : X xy ^ {^li 0, 1, . . . , p— 2} as 
7r(x, y) := max{-l, 0, 1, . . . , p~ 2} subject to j,) n /i^^y) 7^ 

where p is as in Assumption [TJ 

5 The result 

Below is our main theorem. 

Theorems Suppose Assumptions\Tl(^hold. Then system 

s^^^fimli,,y)^h-\y)) (8) 

is a deadbeat observer for system 
Proof. We claim the following 

X e Mti =^ x+e [f{x)]+ (9) 

for all ^ e {0, 1, . . . , p — 1}. Let us prove our claim. Note that x G [x]'^_i yields 



X „ 1 = X 



J by Assumption [2] Since [x]'^_j_ 7^ we have [x]'j'_-^^ [x]o 7^ 
and, consequently, [x]'^_i n [a;]o 7^ 0- Remark[2]then yields [x]^f^^ C 
Starting from we can proceed as 



= fimU,H=^))'^h-\h{x))) 
C f{[x]+_^ n [a;]o) 

= /(HtinHo) 

= /(H^) (10) 
= {f{^)]t- 

Hence ([9]) holds. In particular, (fTO|) gives us 

£ G [0;]+ 1 =^ x+e fi[x]i) (11) 

for all G {0, 1, . . . , p — 1}. Note that x G holds for all x, x. Therefore 

([5]) and Remark [21 imply the existence of ^* G {0, 1, . . . , p — 1} such that 

iPik, X, x) G [<^(fc, x)]+_2 (12) 

for all k > £*. Also, Assumption [1] yields us 

[0(fc, = 0(fc, x) (13) 
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for all fc > p — 1 . Combining ([TT|) , (IT^ , and ([T^ we can write 

X, x) = (/)(fc, x) 

for all fc > p. Hence the result. ■ 

Corollary 1 Consider linear system (g]) with A e M"^" and C G M^'*". ^wp- 
pose pair (C, ^) is observably. Let Sk for k = 0, 1, ... be defined as in 
Lemma{^ Then system 

x+ = A{{x + ASn-2) r\{x + Sq)) 

is a deadbeat observer for system ^ . 



6 System with input 

In this section we look at the case where the evolution of system to be observed 
is dependent not only on the initial condition but also on some exogenous signal, 
which we call the input. To construct a deadbeat observer for such system we 
again make use of sets. 
Consider the system 

x+ = f{x, u) (14a) 
y - h{x) (14b) 

where u Cz U d is the input or some known disturbance (e.g. time). Let 
u = (uo, Ml, . . .), Uk € U, denote an input sequence. The solution of system (|14p 
at time k, starting at initial condition x € X, and having evolved under the 
influence of input sequence u is denoted by (j){k, x, u). Note that (j){0, x, u) — x 
and (p{k + 1, x, u) = f{4>{k, x, u), Uk) for all x, u, and k. 
Now consider the following cascade system 

x^ = f{x, u) (15a) 
x^ G g{x, h{x), u) (15b) 

We denote a solution of subsystem (|15b ) by ip{k, x, x, u). We then have 
•0(0, X, X, u) = X and ip{k + 1, x, x, u) S g{^{k, x, x, u), h{(j){k^ x, u)), Uk) 
for all x, x, u, and k. 

Definition 3 Given g: XxyxU^X, system 

G 9{x, y, u) 

is said to be a deadbeat observer for system (|14p if there exists p > 1 such that 
solutions of system (|15|) satisfy 

ip{k, X, X, u) = (/)(fc, X, u) 

for all X, X, u, and k> p. 

iThat is, rank [C^ A^CT . . . ^(""1)^0^] = n. 
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How to define sets [x]k and [x]^ for system is obvious. We again let 
[x]o ■.= h-\h{x)). 



and (for fc = 0, 1, . . .) 



[a;]fc+i [x]^ n [x\q 



this time with 



[x] 



+ 



U fMk,u). 



The following result is a generalization of Theorem [3] (The demonstration is 
parallel to that of Theorem [3] and hence omitted.) 

Theorem 4 Suppose AssumptionslTl[^ hold. Then system 



is a deadheat observer for system p4p . 

7 Examples 

Here, for two third order nonlinear systems, we construct deadbeat observers. In 
the first example we study a simple autonomous homogeneous system and show 
that the construction yields a homogeneous observer. Hence our method may be 
thought to be somewhat natural in the vague sense that the observer it generates 
inherits certain intrinsic properties of the system. In the second example we aim 
to provide a demonstration on observer construction for a system with input. 

7.1 Homogeneous system 

Consider system ^ with 



where x = [xi X2 x^]"^ . Let X — R^ and 3^ = R. If we let dilation Aa be 



x+ e f{[x] 



+ 

tt(x., y) 



nh^^y), u) 




and h{x) := xi 



A 
A 
A3 



with A G M, then we realize that 



fiAxx)^Axf{x) and h{Axx) ^ Xh{x) . 
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That is, the system is homogeneous [TU] with respect to dilation A. Before 
describing the relevant sets [x]k and [x]~^ we want to mention that / is bijective 
and its inverse is 



Since h{x) — xi we can write 



(a;3 - x?)V3 



Xi 



(16) 





• 


Xi 








a 


: a, /? e laj 






J _ 





(17) 



By ([T6| we can then proceed as 



■(a:3-x?)V3 



/ 



= / 



7 
6 



ix3 - X?)V3 



7 
(5 



7, ^ e R 
: 7, (5 e 



7 

^1/3 

— a;? + 7^ 



7, (5 e 



(18) 



Recall that [a;]i = [x]q n[a;]o. Therefore intersecting sets (fTTj) and (fTSl) we obtain 





• 


Xi 




a;]i = < 




a 


: a € ftj 











We can now construct [x]'^ as 

■(x3-:r3)V3 



/ 



7 
7 



: 7 e 



76 



7 

— + 73 



7€ 



(19) 
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Now note that sets (|17p and HH) intersect at a single point. In particular, 
[x]2 = [x]^ n [x]o — X. Therefore Assumption [1] is satisfied with p = 3. Observe 
also that 



[x]tnh-\y) 



1 






• 


V 




i2 


: 7 e M 


H 




a 


: a, /3 e 


x^ — x\+ 7^ 








.P. 





y 

^3 - ^ 1 + 



which means that ■k{x, y) ^ p — 2 = \ for all x and ?/. The dynamics of the 
deadbeat observer then read 

x+ = f{[x]tnh-\y)) 

{x3-xl + y^y/^ 
xl+y^ 

We finally notice that 

/([A,x]+ n h-\Xy)) - A^f{[x]+ n h-\y)) . 

That is, the deadbeat observer also is homogeneous with respect to dilation A. 

7.2 System with input 

Our second example is again a third order system, this time however with an 
input. Consider system ([T4|) with 



f{x, u) 



X1X2X3 
xz/xi 

^X\X2U 



and h(x) :— x\ . 



Let X = M>Q, 3^ = K>Oj and U — M>o- Let us construct the relevant sets [x]}. 
and [x]^ . We begin with [x]o- 



(20) 





• 


Xl 




a;]o = < 




a 


: a, /3 > oj 






J _ 





Note that / satisfies the following 
/ 



Xiu/{X2X1) 
X2Xy{xiu'^) 



xiu/x\ 



^ U \ — X 
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for all X and u. Hence we can write 



X2xi/{xiu'^) 



= U/ 
= U/ 



ueu 



xij6 
X2S 



X\ulx\ 

Xiu/ {X2X\) 

7 
5 



Xiu/{X2X1) 
X2xl'-f / {xiU^) 

Xiu5/x1 

Xiu/{X2XI) 
X2X\^ / {xiU^) 

Xiu5/x1 



: 7, J > 



7, 5 > 1 , li 
: 7, J > o| , w 
,u \ : 7, (5 > 



Since [x]i — [x\q f^ [xjo, intersecting sets (f20l) and (f2T|) we obtain 





• 






x]i = < 






: a > oj 






2:30; 





We can now construct \x\^ as 

[At = U/ 
- U/ 



a;iu/(x2X3) 
a;2a;3/(a;iu2) 
xiu/xl 



a;iM/(x2X3) 
X2xy {xiu^^"^) 
xiu-f/xl 



7 > > , u 





a;i/7 




! 


X2I 


:7>o| 




.^3/7. 





(21) 



(22) 



Now note that sets (|20p and ([221) intersect at a single point. In particular, 
[2^] 2 — [x\l ^ — X . Therefore Assumption [T] is satisfied with p = 3. Observe 
also that 



[x]+nh-\y) = 







• 


y 




X21 


: 7 > o| n < 




a 


: a, /3 > oj 


2:3/7. 






.P. 





y 

xiX2/y 
xsy/xi 
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which means that 7r(x, y) = p — 2 = 1 for all x and y. The dynamics of the 
deadbeat observer then read 

x+ = f{[x]+nh-\y),u) 
= x-i/xi 

\Jx\XiU 



8 An algorithm for deadbeat gain 

In this section we provide an algorithm to compute the deadbeat observer gain 
for a linear system with scalar output. (The algorithm directly follows from 
Corollary [TJ) Namely, given an observable pair (C, A) with C E M^^" and 
A S M"^", we provide a procedure to compute the gain i e M"^i that renders 
matrix A~ LC nilpotent. Below we let null(-) be some function such that, given 
matrix AI e jj'nx" whose dimension of null space is k, null(Af) is some n x k 
matrix whose columns span the null space of M. 

Algorithm 1 Given C E M.^^'"' and A £ M"^", the following algorithm gener- 
ates deadbeat gain L G M"'*-'^. 

X = null(C) 

for i = 1 : n — 2 

C 



X = null 



null((AX) 



T\T 



end 

AL 



L = 



pre 



CLpYe 

For the interested reader we below give a MATLAB code. Exploiting Algo- 
rithm [1] this code generates a function (which we named dbLfun) whose inputs 
are matrices C and A. The output of the function, as its name indicates, is the 
deadbeat gain L. 

function L = dbLfun(C,A) 

X = null(C) ; 

for i = l:length(A)-2 

X = null([C;null((A*X) ') ']) ; 

end 

Lpre = A*X; 

L = A*Lpre/(C*Lpre) ; 

One can also use the built-in MATLAB function acker to compute the deadbeat 
gain. We can therefore compare dbLfun with acker via a numerical experiment. 
Table[l]gives the experimental results. Number n is the dimension of the system 
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Table 1: Percentages of cases where dbLfun performed better than acker. 



n=3 


n=4 


n=5 


n=6 


n=7 


n=8 


n=9 


n=10 


%51 


%60 


%67 


%74 


%80 


%85 


%87 


%91 



(that is, the number of columns of A matrix) and the numbers at the bottom row 
are the percentages of the cases (among 10^ random trials for each n) in which 
dbLfun performed better than acker. How wc determine which one is better 
in a given case is as follows. Given pair (C, A), we let Li be the gain resulting 
from dbLfun(C,A) and L2 be the gain given by acker (A' ,C' ,zeros(n,l)) 
Then we compare norms |(A — LiC)"| and \ {A — L2C)'^\, neither of which is zero 
due to round-off errors. The function yielding the smaller norm is considered 
to be better. 

9 Conclusion 

For nonlinear systems a method to construct a deadbeat observer is proposed. 
The resultant observer can be considered as a generahzation of the linear dead- 
beat observer. The construction makes use of sets that are generated iteratively. 
Through such iterations, observers are derived for two academic examples. Also, 
for computing the deadbeat gain for a linear system with scalar output, an al- 
gorithm that works no worse than an already existing one is given. 
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